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We search for decays of the type B? s n — > h + h'~ (where h,h' = K or 7r) in 180 pb _1 of pp 
collisions collected at the Tevatron by the upgraded Collider Detector at Fermilab. We report the 
first observation of the new mode -6°— > K + K~ with a yield of 236 ± 32 events, corresponding to 
(fs/fd) x B(B° S -> K+K-)/B(B° -» K + Tv~) = 0.46 ± 0.08(stat.) ± 0.07(syst.), where f s /f d is the 
ratio of production fractions of B® and B°. We find results in agreement with world averages for 
the B° modes, and set the following new limits at 90% CL: B(B°-> K~tx + ) < 5.6 x 10" 6 and 
B{B° S -+ TV + TT-) < 1.7 x 10" 6 



PACS numbers: 13.25.Hw 14.40.Nd 



The decay modes of B mesons into pairs of charm- 
less pseudoscalar mesons are effective probes of the 
quark-mixing (Cabibbo-Kobayashi-Maskawa, CKM) ma- 
trix and are sensitive to potential new physics effects. 
Their branching fractions and CP asymmetries can be 



predicted with good accuracy and compared to rich ex- 
perimental data available for B + and B° mesons, pro- 
duced in large quantities in T(45) decays Q. Measure- 
ments of similar modes predicted, but not yet observed, 
for the B® meson are important to complete our under- 
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standing of B meson decays. The measurement of observ- 
ables from both strange and non-strange B mesons allows 
a cancellation of hadronic uncertainties, thus enhancing 
the precision of the extraction of physics parameters from 
experimental data [1, H, 0, [Bj] . 

The branching fraction of the B®—> K + K~ mode is 
a candidate for observing an unusually large breaking 
of U-spin symmetry, and is sensitive to anomalous elec- 
troweak penguin contributions from new physics 0, S 0] • 
A combination of B ^ ir + ir~ and B®^ K + K~ observ- 
ables has been proposed as a way to directly determine 
the phase of the V u b element of the CKM matrix (angle 
7), or alternatively as a test of our understanding of dy- 
namics of B hadron decays, when compared with other 
determinations of 7 0. The -Bg— > K~ir + mode can also 
be used in measuring 7 [3], and its CP asymmetry is 
a powerful model-independent test [§| of the source of 
the direct CP asymmetry observed in the B ^ K + tt~ 
mode [loj . The B®^ ir + n~ mode proceeds only through 
annihilation diagrams, which are currently poorly known 
and a source of significant uncertainty in many theoret- 
ical calculations [llj. Its features are similar to the as 
yet unobserved B ^ K + K~ mode, but it has a larger 
predicted branching fraction [ll|, [l2| ; a measurement of 
both modes would allow a determination of the strength 
of penguin-annihilation 0] ■ 

In this letter we report the first observation of the de- 
cay B° s — > K + K~ and perform the first measurement in 
hadron collisions of partial widths of B9 •. decays to pairs 
of charged pions and kaons. Throughout this paper, C- 
conjugate modes are implied and branching fractions in- 
dicate CP-averages unless otherwise stated. 

The measurements have been performed in a sample 
of 180 pb _1 of pp collisions at y/s = 1.96 TeV, recorded 
at the Tevatron collider by the upgraded Collider De- 
tector at Fermilab (CDF II). CDF II is a multipurpose 
magnetic spectrometer surrounded by calorimeters and 
muon detectors [HI]. The components of the detector 
pertinent to this analysis are described briefly below. A 
silicon microstrip detector (SVX II) [l4| and a cylindrical 
drift chamber (COT) [3] immersed in a 1.4 T solenoidal 
magnetic field allow reconstruction of charged particles 
in the pseudorapidity range | 77 |< 1.0 [lij ]. The SVX II 
consists of five concentric layers of double-sided silicon 
detectors with radii between 2.5 and 10.6 cm, each pro- 
viding a measurement with 15 /xm resolution in the (f> di- 
rection. The COT has 96 measurement layers, between 
40 and 137 cm in radius, organized into alternating axial 
and ±2° stereo superlayers. The transverse momentum 
resolution is a PT /pr — 0.15%pT/(GeV/c). The specific 
energy loss (dE / dx) of charged particles in the COT can 
be measured from the collected charge, which is encoded 
in the output pulse-width of each wire. 

Data were collected by a three-level trigger system, us- 
ing a set of requirements dedicated to B hadron decays 
into charged particle pairs. At Level 1, charged particle 



tracks are reconstructed in the COT transverse plane by a 
hardware processor (XFT) [17| . Two opposite-curvature 
tracks are required, with reconstructed transverse mo- 
menta pti,PT2 > 2 GeV/c, the scalar sum pxi + Pti > 
5.5 GeV/c, and a transverse opening-angle A</> < 135°. 
At Level 2, the Silicon Vertex Trigger (SVT) [l8] com- 
bines XFT tracks with SVX II hits to measure the im- 
pact parameter d (distance of closest approach to the 
beam line) of each valid track. The requirement of two 
tracks with 100 /mi < d < 1.0 mm reduces light-quark 
background by two orders of magnitude while preserv- 
ing ~ 50% of the signal. A tighter opening angle cut: 
20° < Acf) < 135° selects two-body B decays from multi- 
body with 97% efficiency and reduces background fur- 
ther. Each track pair is then used to form a B candidate, 
which is required to have an impact parameter relative 
to the beam axis ds < 140 /jm and to have travelled a 
transverse distance L xy > 200 /.tm. At Level 3, a farm of 
computers confirms the selection with a full event recon- 
struction. The overall acceptance of the trigger selection 
is ~ 2% for B mesons of px > 4 GeV/c. 

In the offline analysis, combinatoric and light-quark 
backgrounds are effectively rejected by requiring the B 
candidate to be isolated. The isolation cut (I > 0.5) [l9| 
has been chosen, together with tightened cuts on kine- 
matic observables (L xy > 300 (Am, dB < 80 //m, and 
d > 150 (Am), by maximizing the quantity S/(S+B) 1 / 2 
over all possible combinations of cuts. The background 
B is estimated from the data sidebands. The expected 
signal yield S is obtained from a detailed detector simu- 
lation assuming the momentum distribution of B-mesons 
measured by CDF [2fJ, and normalized to the yield ob- 
served after the trigger selection. The overall efficiency 
of the chosen offline selection is ~ 50%. 

No more than one B meson candidate per event sur- 
vives the selection, and a mass is assigned to each, as- 
suming the pion mass for both decay products. The mass 
distribution, shown in Fig. [TJ exhibits an obvious peak in 
the B® s ^ mass region. A binned fit of a Gaussian over an 
exponential background provides an estimate of 893 ± 47 
signal events, with a width a = 38 ± 2 MeV/c 2 , com- 
pared to an expected mass resolution a = 28 MeV/c 2 
for an individual B9 g <. — ► h + h'~ mode. This indicates 
the presence of at least two distinct final states. Siz- 
able signal contributions are expected from two known 
B° modes, B ^ ir + ir~ and B ^ K + ir~ , and two as yet 
unobserved B° s modes, B®— > K + K~ and B®— > K~ir + . 
Figure Q] shows that, as expected, the different modes 
are too closely spaced in mass to be clearly resolved 
and appear instead as a single peak somewhat broader 
than the mass resolution. In addition to mass resolu- 
tion, we use kinematic information along with particle 
identification to extract the different contribution. We 
incorporate all information in an unbinncd likelihood 
fit, to statistically determine the contribution of each 
mode, and the CP asymmetry of the B Q ^ K + tt~ mode 
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A CP = [N(B) - N(B)]/[N(B) + N(B)}. 

For the kinematic portion, we use two loosely corre- 
lated observables to summarize the information carried 
by all possible values of invariant mass of the candidate 
B, resulting from different mass assignments to the two 



outgoing particles [2l|. They are the mass M vv calcu- 



lated with the pion mass assignment to both particles, 
and the signed momentum imbalance a — (1 — Px/pijqi, 
where p\ (j>2) is the lower (higher) of the particle mo- 
menta, and q\ is the sign of the charge of the particle of 
momentum p\. Using these two variables, the mass of 
any particular mode can be expressed, in the relativistic 
limit, as: 



M = M z 

771^7712 7T7T 



(2-\a\)(m 



[i-(M 



(1) 

where mi (7712) is the mass of the lower (higher) momen- 
tum particle (Fig. [2l left). Particle identification (PID) 
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FIG. 1: Invariant mass distribution of B? s s — > h + h'~ candi- 
dates passing all selection requirements, using a pion mass 
assumption for both decay products. Cumulative projections 
of the likelihood fit for each mode are overlaid. 




« = (1 -p/p 2 )xq, 



8 -6 -4 -2 2 4 6 8 10 
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FIG. 2: (left) Average M^n versus a for simulated samples 
of B° events, where K + ir~~ and K~ir + are treated separately. 
The solid curves are the corresponding first-order expressions 
from Eq. fTJ. The corresponding plots for the B® are similar, 
but shifted for the mass difference, (right) Distribution of 
dE/dx (mean COT pulse- width) around the average pion re- 
sponse, for calibration samples of kaons and pions (see text). 



information is provided by the measured dE/dx of the 
two tracks. In order to account for their dependence on 
particle momentum, we include in our fit the scalar sum 
Ptot — Pi + P2 asa fifth observable, which in conjunction 
with a provides unique identification of the momenta of 
both particles. 

With the chosen observables, the likelihood contribu- 
tion of the i th event is written as: 



Ci = (1 



3 



f c 



kin^PID 



1,,-kin /-PID 
oi -bck i -bck 



(2) 



'bck' labels background-related quan- 
j runs over the eight distinguishable 



where the index 
tities, the index 
B9-, — > h + h!~ modes (Fig. [2]), and fj are their fractions, 
to be determined by the fit together with the background 
fraction b. The £^ ln is given by the product of the condi- 
tional probability density of M vn for given a and the joint 
probability distribution Pj(a,ptot)- The mass distribu- 
tion is a Gaussian centered at the value of M M obtained 
from Eq. (Q} by setting the appropriate particle masses 
for each decay mode j. The Gaussian width <j M = 28 ±3 
MeV/c 2 was interpolated from the observed widths of 
other two-body decays (D° — > K~tt + , J/i/j — > 
and T — > /i+^ _ ), and the B° and B° masses are set to 
the values measured by CDF [22| to cancel the common 
systematic uncertainty. The background mass distribu- 
tion is fitted to an exponential function plus a constant. 
The Pj(a,pt t) is parameterized for each mode j by a 
product of polynomial and exponential functions, fitted 
to Monte Carlo samples produced by a detailed detector 
simulation, while the corresponding distribution for the 
background is obtained from the mass sidebands of data. 

The dE/dx response was calibrated over the track- 
ing volume and time by means of a 97%-pure sample 
of 3 x 10 5 D*+ D°7r+ -> [K-TT+]n+ decays, where the 
D° decay products are identified by the charge of the D* + 
pion [23J. The observed response (Fig. [21 right) is well- 
modeled by the convolution of a single-particle response 
function with a common baseline fluctuation, causing 
a 10% correlation between particles in the same event. 
Both effects are quasi- Gaussian with small tails, and have 
been accurately modeled in £ PID . The separation be- 
tween pions and kaons in the range 2 < px < 10 GeV/c 
is nearly constant at 1.4 standard deviations, correspond- 
ing to a resolution 1.7 times worse than a "perfect" PID, 
when measuring the relative fractions of the two parti- 
cles in any given sample. The term allows for in- 
dependent pion, kaon, proton, and electron components, 
which arc free to vary independently in three mass re- 
gions (left, under, and right of the signal peak) to al- 
low for possible variations due to the contribution of 
partially-reconstructed B hadrons in the lower-mass re- 
gion. Muons are indistinguishable from pions with the 
available dE/dx resolution. 

The fit of the data sample returns the yields listed in 
TablelU The observed resolutions are compatible with ex- 
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pectations from fitting Monte Carlo samples of the same 
size. Significant signals are seen for B°— > ir + ir~ , B°^ 
K + ir~ , and the previously unobserved -B°^ K + K~ 
mode, while no evidence is obtained for B®~- ► K~Tt + , 
B®^ 7r + 7r _ , or B°^r K + K~ . As a check of our re- 
sults, we performed an alternative fit based solely on 
kinematical information. Since the B°— > ir + ir~ mode 
is indistinguishable from B®—* K + K~ in absence of PID 
information, we constrain its rate to its world-average 
value [24j ]. This fit confirms the main results, returning 
a yield of 193 ± 55 B®— > K + K~ events. To convert raw 
yields into relative branching fractions, we apply correc- 
tions for the different efficiencies of trigger and offline 
selection requirements for different decay modes; the rel- 
ative efficiency corrections between modes do not exceed 
19%. Most corrections are determined from the detailed 
detector simulation, with the following exceptions which 
are measured using data: A momentum-averaged rela- 
tive isolation efficiency between B° s and B° of 1.07±0.11 
has been determined from fully-reconstructed samples of 
J/tp<p, B° s ^ D~TT + , 5°-* J/ipK*°, and B ^ 
D~n + . The lower specific ionization of kaons with re- 
spect to pions in the COT is responsible for a ~ 5% 
lower efficiency to reconstruct a kaon by the XFT. This 
effect is measured in a sample of D + — > K~-k + tt + decays 
triggered on two tracks, using the unbiased third track. 
The only correction needed by Aqp is a (1.0 ± 0.25)% 
shift due to the different probability for K + and K~ 
to interact with the tracker material. The accuracy of 
our control over instrumental charge asymmetries is con- 
firmed by the smallness of the asymmetry (< 0.5%) mea- 
sured in the D° -> K~n + mode 251. The B"^ K+K~ 



and B 







7r modes require a special treatment, since 



they contain a superposition of the flavor eigenstates of 
the B®. Their time evolution might differ from the flavor- 
specific modes if the width difference Ar s between the 
B® mass eigenstates is significant. The current result is 
derived under the assumption that both modes are dom- 
inated by the short-lived B® component , that T s = T^, 
and Ar s /r s = 0.12 ± 0.06 j2f|. The latter uncertainty is 
included in estimating the overall systematic uncertainty. 

The dominant contributions to the systematic uncer- 
tainty are: the statistical uncertainty on isolation effi- 
ciency (B® modes) , possible charge asymmetry of back- 
ground (Aqp), and final state photon radiation > 
7r + 7r~). The latter is conservatively estimated with the 
full effect predicted by QED calculations [St]]. Smaller 
systematic uncertainties are assigned for: mass scale and 
resolution; dE/dx response model ; trigger efficiencies; 
background shape and kinematics; B meson masses, life- 
times, and differences in momenta, allowed to vary by 
a factor (mgo — m#o)/m#o due to fragmentation ef- 
fects (28| . 

The relative branching fractions obtained after apply- 
ing all corrections are listed in Table HI where fd and 
f a indicate the production fractions respectively of B° 



and B® from fragmentation of a & quark in pp collisions. 
Upper limits are quoted for modes in which no signifi- 
cant signal is observed [23]. We also list absolute results 
obtained by normalizing our data to the world-average 
of B(B ^ K + ir~) and assuming for f s /fd the world- 



average from pp and e + e experiments 24 



The rate of the newly observed mode B®—* K + K~ 
favors the higher value (36 ± 7) x 10 -6 predicted by cal- 
culations based on QCD sum rules [3, E| implying large 
U-spin breaking in this process, although it is not statisti- 
cally incompatible with the expectation B{B®^ K + K~) 
= B(B ^ K + ir~) from the assumption of exact U- 
spin symmetry and negligible spectator contributions. 
We also derive the ratio of U-spin-conjugate decays: 
(f d /f s ) x B(B° - 7r+7r-)/B(B Q s - K+K~) = 0.45 ± 
0.13 ± 0.06, which can be related to the CP asymmetries 
in the tt+tt" mode and to the CKM angle 7 @]. 

Our results for the B° are in agreement with world av- 
erage values: B(B ^ Tr+Tr") = (4.6 ± 0.4) x 10~ 6 and 
A C p(B a ^ K+tt-) = -0.113 ± 0.020 0] although our 
Acp measurement is also compatible with zero. The limit 
set on B®— > K~ir + indicates a value at the lower end 
of current expectations [!, HH . The limit for the anni- 
hilation mode 7T + 7r - is a large improvement over 
the previous best limit [31|. approaching the expectations 
from recent calculations [121, K32j • 



In summary, we have measured relative branching frac- 
tions of B? s s mesons into pairs of charmless charged 
mesons. We find results in agreement with current world 
averages for B° modes and observe for the first time the 
Bg— » K + K~ mode. We set upper limits on unobserved 
modes B°^ K+K~, B°^ R-1T+, and 5°-> ir+ir-. 
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TABLE I: Summary of results. The yields of the two annihilation modes (last two rows) were fixed to zero when fitting 
for the four main modes. Absolute branching fractions are normalized to the the world-average values B(B°— > K + iv~) = 
(18.9 ± 0.7) x 10 -6 and fs/fd = 0.26 ± 0.039 [24| . The first quoted uncertainty is statistical, the second is systematic. 



Mode Yield 


Measured Quantity 


Derived B {10~") 


B u -> K+n~ 542 ± 30 

B°^TT + TT~ 121 ± 27 

B ( > -> K+R- 236 ± 32 
B°^K-tt+ 3 ± 25 


A C p = -0.013 ± 0.078 ± 0.012 

ggr^+r-) = °- 21 ± °-° 5 ± °- 03 

b*$fc£+V-] = 0-46 ± 0.08 ± 0.07 
^ITJCkIC] <0.08@90%CL 


3.9 ± 1.0 ± 0.6 

33 ± 6 ± 7 
< 5.6 @ 90% CL 


B°->tt + tt- -10 ± 15 
B° -> K+R" 10 ± 23 


<0.05@90%CL 
Itfel^ < 0.10 O 90% CL 


< 1.7 @ 90% CL 

< 1.8 @ 90% CL 
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